Exome sequencing followed by large-scale genotyping suggests a limited role for moderately rare risk factors of strong effect in schizophrenia.

Schizophrenia is a severe psychiatric disorder with strong heritability and marked heterogeneity in symptoms, course, and treatment response. There is strong interest in identifying genetic risk factors that can help to elucidate the pathophysiology and that might result in the development of improved treatments. Linkage and genome-wide association studies (GWASs) suggest that the genetic basis of schizophrenia is heterogeneous. However, it remains unclear whether the underlying genetic variants are mostly moderately rare and can be identified by the genotyping of variants observed in sequenced cases in large follow-up cohorts or whether they will typically be much rarer and therefore more effectively identified by gene-based methods that seek to combine candidate variants. Here, we consider 166 persons who have schizophrenia or schizoaffective disorder and who have had either their genomes or their exomes sequenced to high coverage. From these data, we selected 5,155 variants that were further evaluated in an independent cohort of 2,617 cases and 1,800 controls. No single variant showed a study-wide significant association in the initial or follow-up cohorts. However, we identified a number of case-specific variants, some of which might be real risk factors for schizophrenia, and these can be readily interrogated in other data sets. Our results indicate that schizophrenia risk is unlikely to be predominantly influenced by variants just outside the range detectable by GWASs. Rather, multiple rarer genetic variants must contribute substantially to the predisposition to schizophrenia, suggesting that both very large sample sizes and gene-based association tests will be required for securely identifying genetic risk factors. © 2012 The American Society of Human Genetics

Screening the human exome: a comparison of whole genome and whole transcriptome sequencing

BACKGROUND: There is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important. RESULTS: Here we provide a systematic exploration of how well RNA-Seq can identify human coding variants by comparing variants identified through high coverage whole-genome sequencing to those identified by high coverage RNA-Seq in the same individual. This comparison allowed us to directly evaluate the sensitivity and specificity of RNA-Seq in identifying coding variants, and to evaluate how key parameters such as the degree of coverage and the expression levels of genes interact to influence performance. We find that although only 40% of exonic variants identified by whole genome sequencing were captured using RNA-Seq; this number rose to 81% when concentrating on genes known to be well-expressed in the source tissue. We also find that a high false positive rate can be problematic when working with RNA-Seq data, especially at higher levels of coverage. CONCLUSIONS: We conclude that as long as a tissue relevant to the trait under study is available and suitable quality control screens are implemented, RNA-Seq is a fast and inexpensive alternative approach for finding coding variants in genes with sufficiently high expression levels

Streptomyces‐derived quorum‐sensing systems engineered for adjustable transgene expression in mammalian cells and mice

Prokaryotic transcriptional regulatory elements have been adopted for controlled expression of cloned genes in mammalian cells and animals, the cornerstone for gene‐function correlations, drug discovery, biopharmaceutical manufacturing as well as advanced gene therapy and tissue engineering. Many prokaryotes have evolved specific molecular communication systems known as quorum‐sensing to coordinate population‐wide responses to physiological and/or physicochemical signals. A generic bacterial quorum‐sensing system is based on a diffusible signal molecule that prevents binding of a repressor to corresponding operator sites thus resulting in derepression of a target regulon. In Streptomyces, a family of butyrolactones and their corresponding receptor proteins, serve as quorum‐sensing systems that control morphological development and antibiotic biosynthesis. Fusion of the Streptomyces coelicolor quorum‐sensing receptor (ScbR) to a eukaryotic transactivation domain (VP16) created a mammalian transactivator (SCA) which binds and adjusts transcription from chimeric promoters containing an SCA‐specific operator module (PSPA). Expression of erythropoietin or the human secreted alkaline phosphatase (SEAP) by this quorum‐sensor‐regulated gene expression system (QuoRex) could be fine‐tuned by non‐toxic butyrolactones in a variety of mammalian cells including human primary and mouse embryonic stem cells. Following intraperitoneal implantation of microencapsulated Chinese hamster ovary cells transgenic for QuoRex‐controlled SEAP expression into mice, the serum levels of this model glycoprotein could be adjusted to desired concentrations using different butyrolactone dosing regime

Prioritization of Epilepsy Associated Candidate Genes by Convergent Analysis

Epilepsy is a severe neurological disorder affecting a large number of individuals, yet the underlying genetic risk factors for epilepsy remain unclear. Recent studies have revealed several recurrent copy number variations (CNVs) that are more likely to be associated with epilepsy. The responsible gene(s) within these regions have yet to be definitively linked to the disorder, and the implications of their interactions are not fully understood. Identification of these genes may contribute to a better pathological understanding of epilepsy, and serve to implicate novel therapeutic targets for further research.In this study, we examined genes within heterozygous deletion regions identified in a recent large-scale study, encompassing a diverse spectrum of epileptic syndromes. By integrating additional protein-protein interaction data, we constructed subnetworks for these CNV-region genes and also those previously studied for epilepsy. We observed 20 genes common to both networks, primarily concentrated within a small molecular network populated by GABA receptor, BDNF/MAPK signaling, and estrogen receptor genes. From among the hundreds of genes in the initial networks, these were designated by convergent evidence for their likely association with epilepsy. Importantly, the identified molecular network was found to contain complex interrelationships, providing further insight into epilepsy's underlying pathology. We further performed pathway enrichment and crosstalk analysis and revealed a functional map which indicates the significant enrichment of closely related neurological, immune, and kinase regulatory pathways.The convergent framework we proposed here provides a unique and powerful approach to screening and identifying promising disease genes out of typically hundreds to thousands of genes in disease-related CNV-regions. Our network and pathway analysis provides important implications for the underlying molecular mechanisms for epilepsy. The strategy can be applied for the study of other complex diseases

Detection of Pesticides in Active and Depopulated Beehives in Uruguay

The influence of insecticides commonly used for agricultural purposes on beehive depopulation in Uruguay was investigated. Honeycombs, bees, honey and propolis from depopulated hives were analyzed for pesticide residues, whereas from active beehives only honey and propolis were evaluated. A total of 37 samples were analyzed, representing 14,800 beehives. In depopulated beehives only imidacloprid and fipronil were detected and in active beehives endosulfan, coumaphos, cypermethrin, ethion and chlorpyrifos were found. Coumaphos was present in the highest concentrations, around 1,000 μg/kg, in all the propolis samples from active beehives. Regarding depopulated beehives, the mean levels of imidacloprid found in honeycomb (377 μg/kg, Standard Deviation: 118) and propolis (60 μg/kg, Standard Deviation: 57) are higher than those described to produce bee disorientation and fipronil levels detected in bees (150 and 170 μg/kg) are toxic per se. The other insecticides found can affect the global fitness of the bees causing weakness and a decrease in their overall productivity. These preliminary results suggest that bees exposed to pesticides or its residues can lead them in different ways to the beehive

Alternating hemiplegia of childhood: evolution over time and mouse model corroboration

Alternating hemiplegia of childhood is a rare neurodevelopmental disorder caused by ATP1A3 mutations. Some evidence for disease progression exists, but there are few systematic analyses. Here, we evaluate alternating hemiplegia of childhood progression in humans and in the D801N knock-in alternating hemiplegia of childhood mouse, Mashlool, model. This study performed an ambidirectional (prospective and retrospective data) analysis of an alternating hemiplegia of childhood patient cohort (n = 42, age 10.24 ± 1.48 years) seen at one US centre. To investigate potential disease progression, we used linear mixed effects models incorporating early and subsequent visits, and Wilcoxon Signed Rank test comparing first and last visits. Potential early-life clinical predictors were determined via multivariable regression. We also compared EEG background at first encounter and at last follow-up. We then performed a retrospective confirmation study on a multicentre cohort of alternating hemiplegia of childhood patients from France (n = 52). To investigate disease progression in the Mashlool mouse, we performed behavioural testing on a cohort of Mashlool- mice at prepubescent and adult ages (n = 11). Results: US patients, over time, demonstrated mild worsening of non-paroxysmal disability index scores, but not of paroxysmal disability index scores. Increasing age was a predictor of worse scores: P < 0.0001 for the non-paroxysmal disability index, intellectual disability scale and gross motor scores. Earliest non-paroxysmal disability index score was a predictor of last visit non-paroxysmal disability index score (P = 0.022), and earliest intellectual disability score was a predictor of last intellectual disability score (P = 0.035). More patients with EEG background slowing were noted at last follow-up as compared to initial (P = 0.015). Similar worsening of disease with age was also noted in the French cohort: age was a significant predictor of non-paroxysmal disability index score (P = 0.001) and first and last non-paroxysmal disability index score scores significantly differed (P = 0.002). In animal studies, adult Mashlool mice had, as compared to younger Mashlool mice, (i) worse balance beam performance; (ii) wider base of support; (iii) higher severity of seizures and resultant mortality; and (iv) no increased predisposition to hemiplegic or dystonic spells. In conclusion, (i) non-paroxysmal alternating hemiplegia of childhood manifestations show, on average over time, progression associated with severity of early-life non-paroxysmal disability and age. (ii) Progression also occurs in Mashlool mice, confirming that ATP1A3 disease can lead to age-related worsening. (iii) Clinical findings provide a basis for counselling patients and for designing therapeutic trials. Animal findings confirm a mouse model for investigation of underlying mechanisms of disease progression, and are also consistent with known mechanisms of ATP1A3-related neurodegeneration

A Genome-Wide Investigation of SNPs and CNVs in Schizophrenia

We report a genome-wide assessment of single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) in schizophrenia. We investigated SNPs using 871 patients and 863 controls, following up the top hits in four independent cohorts comprising 1,460 patients and 12,995 controls, all of European origin. We found no genome-wide significant associations, nor could we provide support for any previously reported candidate gene or genome-wide associations. We went on to examine CNVs using a subset of 1,013 cases and 1,084 controls of European ancestry, and a further set of 60 cases and 64 controls of African ancestry. We found that eight cases and zero controls carried deletions greater than 2 Mb, of which two, at 8p22 and 16p13.11-p12.4, are newly reported here. A further evaluation of 1,378 controls identified no deletions greater than 2 Mb, suggesting a high prior probability of disease involvement when such deletions are observed in cases. We also provide further evidence for some smaller, previously reported, schizophrenia-associated CNVs, such as those in NRXN1 and APBA2. We could not provide strong support for the hypothesis that schizophrenia patients have a significantly greater “load” of large (>100 kb), rare CNVs, nor could we find common CNVs that associate with schizophrenia. Finally, we did not provide support for the suggestion that schizophrenia-associated CNVs may preferentially disrupt genes in neurodevelopmental pathways. Collectively, these analyses provide the first integrated study of SNPs and CNVs in schizophrenia and support the emerging view that rare deleterious variants may be more important in schizophrenia predisposition than common polymorphisms. While our analyses do not suggest that implicated CNVs impinge on particular key pathways, we do support the contribution of specific genomic regions in schizophrenia, presumably due to recurrent mutation. On balance, these data suggest that very few schizophrenia patients share identical genomic causation, potentially complicating efforts to personalize treatment regimens

Performance of Genotype Imputation for Rare Variants Identified in Exons and Flanking Regions of Genes

Genotype imputation has the potential to assess human genetic variation at a lower cost than assaying the variants using laboratory techniques. The performance of imputation for rare variants has not been comprehensively studied. We utilized 8865 human samples with high depth resequencing data for the exons and flanking regions of 202 genes and Genome-Wide Association Study (GWAS) data to characterize the performance of genotype imputation for rare variants. We evaluated reference sets ranging from 100 to 3713 subjects for imputing into samples typed for the Affymetrix (500K and 6.0) and Illumina 550K GWAS panels. The proportion of variants that could be well imputed (true r2>0.7) with a reference panel of 3713 individuals was: 31% (Illumina 550K) or 25% (Affymetrix 500K) with MAF (Minor Allele Frequency) less than or equal 0.001, 48% or 35% with 0.001<MAF< = 0.005, 54% or 38% with 0.005<MAF< = 0.01, 78% or 57% with 0.01<MAF< = 0.05, and 97% or 86% with MAF>0.05. The performance for common SNPs (MAF>0.05) within exons and flanking regions is comparable to imputation of more uniformly distributed SNPs. The performance for rare SNPs (0.01<MAF< = 0.05) was much more dependent on the GWAS panel and the number of reference samples. These results suggest routine use of genotype imputation for extending the assessment of common variants identified in humans via targeted exon resequencing into additional samples with GWAS data, but imputation of very rare variants (MAF< = 0.005) will require reference panels with thousands of subjects

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition

Contribution of Somatic Ras/Raf/Mitogen-Activated Protein Kinase Variants in the Hippocampus in Drug-Resistant Mesial Temporal Lobe Epilepsy

Importance: Mesial temporal lobe epilepsy (MTLE) is the most common focal epilepsy subtype and is often refractory to antiseizure medications. While most patients with MTLE do not have pathogenic germline genetic variants, the contribution of postzygotic (ie, somatic) variants in the brain is unknown. Objective: To test the association between pathogenic somatic variants in the hippocampus and MTLE. Design, Setting, and Participants: This case-control genetic association study analyzed the DNA derived from hippocampal tissue of neurosurgically treated patients with MTLE and age-matched and sex-matched neurotypical controls. Participants treated at level 4 epilepsy centers were enrolled from 1988 through 2019, and clinical data were collected retrospectively. Whole-exome and gene-panel sequencing (each genomic region sequenced more than 500 times on average) were used to identify candidate pathogenic somatic variants. A subset of novel variants was functionally evaluated using cellular and molecular assays. Patients with nonlesional and lesional (mesial temporal sclerosis, focal cortical dysplasia, and low-grade epilepsy-associated tumors) drug-resistant MTLE who underwent anterior medial temporal lobectomy were eligible. All patients with available frozen tissue and appropriate consents were included. Control brain tissue was obtained from neurotypical donors at brain banks. Data were analyzed from June 2020 to August 2022. Exposures: Drug-resistant MTLE. Main Outcomes and Measures: Presence and abundance of pathogenic somatic variants in the hippocampus vs the unaffected temporal neocortex. Results: Of 105 included patients with MTLE, 53 (50.5%) were female, and the median (IQR) age was 32 (26-44) years; of 30 neurotypical controls, 11 (36.7%) were female, and the median (IQR) age was 37 (18-53) years. Eleven pathogenic somatic variants enriched in the hippocampus relative to the unaffected temporal neocortex (median [IQR] variant allele frequency, 1.92 [1.5-2.7] vs 0.3 [0-0.9]; P =.01) were detected in patients with MTLE but not in controls. Ten of these variants were in PTPN11, SOS1, KRAS, BRAF, and NF1, all predicted to constitutively activate Ras/Raf/mitogen-activated protein kinase (MAPK) signaling. Immunohistochemical studies of variant-positive hippocampal tissue demonstrated increased Erk1/2 phosphorylation, indicative of Ras/Raf/MAPK activation, predominantly in glial cells. Molecular assays showed abnormal liquid-liquid phase separation for the PTPN11 variants as a possible dominant gain-of-function mechanism. Conclusions and Relevance: Hippocampal somatic variants, particularly those activating Ras/Raf/MAPK signaling, may contribute to the pathogenesis of sporadic, drug-resistant MTLE. These findings may provide a novel genetic mechanism and highlight new therapeutic targets for this common indication for epilepsy surgery